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Two-component signal transduction (TCST) is the 
predominant signaling scheme used in bacteria to 
sense and respond to environmental changes in 
order to survive and thrive. A typical TCST system 
consists of a sensor histidine kinase to detect 
external signals and an effector response regula- 
tor to respond to external changes. In the signal- 
ing scheme, the histidine kinase phosphorylates 
and activates the response regulator, which func- 
tions as a transcription factor to modulate gene 
expression. One promising strategy toward anti- 
bacterial development is to target TCST regulatory 
systems, specifically the response regulators to 
disrupt the expression of genes important for viru- 
lence. In Salmonella enterlca, the PhoQ/PhoP sig- 
nal transduction system is used to sense and 
respond to low magnesium levels and regulates 
the expression for over 40 genes necessary for 
growth under these conditions, and more interest- 
ingly, genes that are important for virulence. In 
this study, a hybrid approach coupling computa- 
tional and experimental methods was applied to 
identify drug-like compounds to target the PhoP 
response regulator. A computational approach of 
structure-based virtual screening combined with a 
series of biochemical and biophysical assays was 
used to test the predictability of the computa- 
tional strategy and to characterize the mode of 
action of the compounds. Eight compounds from 
virtual screening inhibit the formation of the 
PhoP-DNA complex necessary for virulence gene 



regulation. This investigation served as an initial 
case study for targeting TCST response regulators 
to modulate the gene expression of a signal 
transduction pathway important for bacterial viru- 
lence. With the increasing resistance of patho- 
genic bacteria to current antibiotics, targeting 
TCST response regulators that control virulence is 
a viable strategy for the development of antimi- 
crobial therapeutics with novel modes of action. 
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Infectious diseases have evolved resistance to most clinical antibi- 
otics. Antibiotic resistance occurs at low levels in natural popula- 
tions, but can become prevalent within a few years of the clinical 
adoption of an antibiotic. Antimicrobial therapeutics currently in 
clinical use either inhibit bacterial growth or induce death. One 
promising strategy is to combat virulence per se without inhibiting 
growth or inducing death, so less selective pressure will cause the 
bacteria to generate resistance. With the emergence of bacterial 
strains resistant to multiple antibiotics, there is an urgent need for 
the development of antibiotics with different modes of action less 
subjective to the development of resistance. 

Two-component signal transduction (TCST) is the predominant sig- 
naling scheme in bacteria to sense and respond to environmental 
changes for survival and proliferation (1-5). TCST regulatory sys- 
tems are modular in terms of their arrangement of domains within 
their proteins within various pathways. In general, TCST regulatory 
systems are comprised of a transmembrane sensor histidine kinase 
(HK) and an intracellular receiver response regulator (RR) with con- 
served sequence, structural, and biochemical properties, allowing 
them to readily adapt to various modes of intracellular signaling. 
These signaling systems typically couple environmental stimuli to 
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an adaptive response, participating in fundamental processes such 
as regulating metabolism, as well as more specialized functions 
such as controlling virulence for the pathogen's host. 

The PhoQ/PhoP two-component regulator system is a major regula- 
tor of virulence in Salmonella enterica serovar Typhimurium and 
also in a number of other gram-negative bacterial pathogens (e.g.. 
Shigella flexneri, Yersinia pestis, Neisseria meningitidis) (6-11). 
PhoQ/PhoP in S. enterica is activated by low extracellular Mg 2+ lev- 
els, acidic pH, and antimicrobial peptides (typical of human gut con- 
ditions during infection) to control various physiological and 
virulence functions (7,12-14). In the signaling cascade (Figure 1), 
the PhoQ histidine kinase is activated by low extracellular magne- 
sium levels and is autophosphorylated at a histidine residue. PhoQ 
subsequently transfers the phosphate group from the conserved his- 
tidine of PhoQ to the conserved aspartate on the PhoP response 
regulator. Phosphorylation of PhoP presumably induces a conforma- 
tional change to mediate homodimerization for DNA binding. The 
PhoP homodimer functions as a transcription factor by recognizing 
and binding to phoP boxes in promoters of PhoP-regulated genes. 
Through this mechanism, PhoP regulates expression of approxi- 
mately 3% of the Salmonella genes in response to low magnesium 
levels to control physiological and virulence functions. The S. enteri- 
ca PhoQ/PhoP signaling pathway is one of the better characterized 
bacterial TCST systems demonstrated to be important for virulence 
regulation. 

PhoP is a member of the OmpR/PhoB structural family, the largest 
structural family making up approximately 30% of all TCST response 
regulators. OmpR/PhoB family members are typically composed of 
two domains connected by a flexible linker: an N-terminal receiver 
domain and a C-terminal DNA-binding domain. Response regulators 
of the OmpR/PhoB family are characterized by a conserved a4-/?5- 
a5 structural motif of the receiver domain, a plastic surface impor- 
tant for dimerization. Owing to the conservation and importance of 
the 0:4-^5-0:5 structural motif for function among OmpR/PhoB family 
members, it was hypothesized that bacterial virulence could be inhib- 



ited by disruption of the PhoQ/PhoP signaling pathway, specifically 
by targeting the s:4-^5-o:5 interface of the PhoP response regulator. 
Preventing this essential protein-protein interaction (PPI) would inhi- 
bit formation of the PhoP-DNA complex and its function as a tran- 
scription factor to regulate gene expression. PhoP was chosen in this 
study as a prototype of the OmpR/PhoB family of response regula- 
tors to probe for PPI 'hot-spots' at the oc4-/J5-a5 interface (Figure 2), 
in efforts to identify and characterize potential binding sites. The 
response regulator N-terminal o<4-/?5-o:5 interface shares a common 
set of hydrophobic and charged residues involved in van der Waals 
contact and salt-bridges important for homodimerization and function 
and is conserved among other response regulators across different 
bacteria species (2,15). High-resolution X-ray crystal structures of 
both inactivated (Protein Data Bank (PDB) ID: 2PKX) and activated 
(PDB ID: 2PL1) Escherichia coli PhoP were available with an interface 
highly similar in sequence to the one in S. enterica (differing only by 
one residue at the o:4-^5-o:5 interface) (16). For these reasons, 
S. enterica PhoP was an attractive target for investigation via struc- 
ture-based drug design to test the effects of response regulator inhi- 
bition and its potential for virulence regulation. 

Targeting bacterial signal transduction systems has only recently 
been demonstrated to be an effective potential strategy for antibi- 
otics development. Rasko et al. (17) targeted the sensor TCST com- 
ponent, the QseC histidine kinase, by the prevention of 
autophosphorylation, which led to disruption of the signaling cas- 
cade important for virulence regulation. Shaknovich ef al. demon- 
strated the feasibility of a small molecule for homodimer inhibition 
and virulence gene regulation when they discovered virstatin to tar- 
get the Vibrio cholerae ToxT (18,19). These studies demonstrated 
the feasibility of drug-like molecules targeting gene expression 
important for virulence regulation as a potential strategy for antibi- 
otics development. 

A prototype of the predominant class of bacterial signal transduc- 
tion important for bacterial virulence is investigated as a proof-of- 
concept study toward this new strategy for antibiotics development. 




Figure 1. Schematic diagram of the Salmonella enterica PhoQ/PhoP two-component signal transduction system. PhoQ is the sensor histi- 
dine kinase, and PhoP is the response regulator. PhoQ senses low extracellular magnesium levels, leading to autophosphorylation at a con- 
served histidine residue. PhoQ transfers to phosphate group to a conserved aspartate residue on PhoP. Phosphorylation of PhoP mediates 
activation by causing a conformational change, allowing PhoP to homodimerize. PhoP recognizes phoP boxes at its DNA promoters (e.g., phoP, 
phoQ, mgtA) and function as a transcription factor to regulate virulence gene expression (including positive feedback loop). 
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Figure 2: Critical salt-bridges at the PhoP a4-/?5-a5 interface 
important for PhoP homodimerization and function. Residues impor- 
tant for dimerization (site-directed mutagenesis studies, Stock, Gao 
& Stock unpublished) are shown in capped sticks and labeled by 
their residue name and number. Mutation of one of these residues 
decreases its ability to homodimerize. The PhoP homodimer (PDB 
ID: 2PKX) (cartoon) binds to phoP boxes in promoters of PhoP-regu- 
lated genes to modulate virulence gene expression. 

TCST systems predominate in control of bacterial expression and 
are completely absent in humans, making them an attractive class 
of targets for the development of new antibiotics with novel modes 



of action. To our knowledge, there are currently no known inhibitors 
of TCST response regulators. Drug-like compounds targeting PhoP, 
specifically the functionally important cz4-/i5-a5 interface, should 
selectively disrupt its function as a transcription factor and inhibit 
the expression of critical virulence genes. In this study, a hybrid 
approach coupling computational and experimental methods 
(Figure 3) was used to predict, validate, and characterize drug-like 
inhibitors of the S. enterica PhoP response regulator. 

Methods and Materials 

Overview of computational strategy 

Figure 3A illustrates an overview of the computational strategy to 
prioritize compounds targeting the ocA-fib-xb interface of PhoP. 
Structure-based virtual screening was used to screen a drug-like 
version of the National Cancer Institute (NCI) Diversity library 
(n = 1420) (20) with the crystal structure of the activated PhoP (PDB 
ID: 2PL1) (16). The docking procedure (labeled 'Docking') was per- 
formed using AutoDock to predict binding poses. A consensus scor- 
ing method (labeled 'Scoring') consisting of the predicted binding 
affinities from AutoDock (21), X-Score (22), and CScore 3 was used 
to better assess the binding affinities of the resulting poses. The 
top-ranked compounds were tested experimentally to assess biolog- 
ical activity and characterize the mode of action. As an alternative 
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Figure 3: Schematic diagram of the computational (A) and experimental (B) workflow to predict and test for PhoP response regulator 
inhibitors. In the computational method (A), a drug-like version of the NCI Diversity library [n = 1420) was screened for inhibitors of activated 
PhoP (PDB ID: 2PL1). Docking was performed using AutoDock 4.0.1 and FRED (OpenEye). The predicted binding poses were rescored using X- 
Score and CScore. A consensus scoring scheme composed of the predict affinities from AutoDock, X-Score, and CScore was used to rank the 
compounds. The top-ranked compounds were tested experimentally. Similarity search (Tanimoto Index) was performed on eight compounds 
with inhibition activity to search for structural analogs. In the experimental strategy (B), the set of prioritized compounds were initially tested 
using EMSA. Compounds displaying inhibition activity by EMSA were further characterized for effects on dimerization using native PAGE, SEC, 
and FRET. Eight compounds inhibited formation of the protein-DIMA complex necessary for virulence gene regulation. 



Chem Biol Drug Des 2012; 79: 1007-1017 



1009 



Tang et al. 



computational method to identify potential PhoP inhibitors, FRED b 
was used for both docking and scoring (consensus method using 
PLP, ChemScore, OEChemScore). The top-ranked compounds from 
FRED [n = 40) were also tested experimentally. 

From the experimental results, a similarity search based on the Tan- 
imoto index (labeled 'Similarity Search') was performed to identify 
structural analogs [n = 36). These analogs were also experimentally 
tested. A total of 255 compounds were tested experimentally and 
verified eight compounds that inhibited formation of the PhoP-DNA 
complex necessary for gene regulation. 



Compound library and PhoP X-ray crystal 
structures 

A drug-like version of the National Cancer Institute (NCI) Diversity I 
library composed of 1420 compounds [derived from a larger version 
of the NCI library of 140 000 compounds filtered based on criteria 
characteristic of drug-like compounds derived from Lipinski's Rule of 
5 (23)] was screened (20,24). Compounds in the NCI Diversity library 
have purity 90% or better as characterized by liquid chromatogra- 
phy-mass spectrometry (LC-MS). The compound library was down- 
loaded in SD format (http://gfscore.cnrs-mrs.fr/download/ 
diversity1440.sdf) and converted to mol2 format using Open Babel 
(http://www.openbabel.org). Inactivated and activated £ coli PhoP 
(PDB ID: 2PKX and 2PL1, respectively) were downloaded from the 
Protein Data Bank (16). 



Structure-based virtual screening 

AUTODOCK 

AutoDock 4.0.1 was used to predict binding poses for compounds in 
the NCI Diversity library for experimental testing (21,25). AutoDock 
was used to convert the ligand structures from the mol2 format to 
the AutoDock pdbqt format, with explicit hydrogen atoms and calcu- 
lated Gasteiger charges (26). AutoDock Tools (ADT) was used to 
prepare the protein structures (27). Polar hydrogens were added to 
the protein target, and Gasteiger charges were assigned. The struc- 
ture files were saved in the pdbqt format. 

Docking was performed to identify low-energy conformations (bind- 
ing poses) of the compounds to sterically and chemically comple- 
ment the binding site. The protein search area (grid spacing of 
0.375 A with dimensions of 40 x 40 x 40 A) was centered at the 
a4-/i5-a5 motif where residues critical for dimerization in PhoP 
(e.g., Arg 111 and Arg 118 of activated PhoP) and other OmpR/ 
PhoB family members are located, as shown in Figure 4 (28). 
Lamarckian genetic algorithm was used to perform the ligand con- 
formational searches to result in 30 binding poses for each ligand. 
The default parameters were used; to increase the sampling and 
accuracy, the following parameters were modified: ga_pop_ 
size = 200; ga_num_evals = 5 000 000; ga_run = 30. 

A drug-like version of the NCI Diversity library (1420 compounds) 
(20) filtered based on drug-like and lead-like features was screened 
using AutoDock 4.0.1 (29). Virtual screening was performed on the 
£ co//-activated PhoP (PDB ID: 2PL1), because the sequence at the 




Figure 4: To prioritize compounds in NCI Diversity (n = 1420) to 
identify PhoP dimerization inhibitors, docking searches were con- 
strained at the interface regions where salt-bridges critical for 
dimerization are located (AutoDock grid boxes used illustrated by 
hashed lines). Residues critical for dimerization (R111 and R118) are 
labeled and shown in CPK, in addition to PhoP (PDB ID: 2PKX) chain 
A (ribbon) and chain B (surface). 



dimer interface was highly similar to the S. enterica PhoP. [S. enter- 
ics PhoP contains SER 93, while £ coli PhoP has ALA 93. Modeling 
was not performed to change the ALA 93 to SER 93 because this 
residue was not present in the specified docking regions.) AutoDock 
was used to generate docked poses of the NCI Diversity compounds 
at the PhoP dimer interface. 



FRED (OpenEye) 

fred was used as an alternative computational tool for both docking 
and scoring. 0MEGA2 was used to generate bioactive ligand con- 
formations according to the parameters from the study by Bostrom 
et al. (30). Default parameters were used for fred docking searches. 
Binding site was marked by using a docked ligand from the NCI 
Diversity library and specifying the regions within 5 A of it. The 
default fred scoring function was used to predict binding affinities 
to rank ligands for experimental testing. 



Rescoring of predicted poses 

Owing to the known limitations of scoring functions to accurately 
predict binding affinities (31-33), especially in the case where the 
compounds are docked to a novel target with a relatively flat bind- 
ing site typical of protein dimer interfaces, a consensus scoring 
scheme was used to prioritize and identify tight-binding ligands 
(31,33-37). Consensus scoring has proven to be effective in increas- 
ing the enrichment rate in a number of studies (37,38). Ligand bind- 
ing poses predicted by AutoDock were rescored using CScore in 
SYBYL 7.3, which is composed of four separate scoring functions 
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(D-score, ChemScore, PMF-score, and G-score), and X-Score (22,39). 
The consensus score for CScore resulted in integer scores between 
1 and 4, with four being the highest. Docked poses were then 
sorted based on the consensus CScore value and subsequently by 
the ChemScore binding affinities. The consensus score estimated by 
X-Score with default parameters, composed of three different scor- 
ing functions (HPScore, HMScore, and HSScore), was used. The top 
15% of the docked poses ranked by each of the scoring functions 
(the AutoDock binding free energy, the CScore consensus method, 
and the X-Score consensus method) were selected, and the docked 
poses found in all three scoring methods (6390 poses representing 
179 compounds) were selected for experimental testing. 

This computational protocol of docking and consensus scoring has 
been demonstrated to be effective in our research group for identi- 
fying lead compounds from the NCI Diversity library for various PPI 
studies (40). Because there are no known inhibitors of PhoP or any 
TCST response regulators with which to validate a computational 
protocol (there are only a small number of structural examples of 
protein-protein interaction inhibitors available as reviewed by Wells 
(41) and no examples to our knowledge with a binding site that is 
structurally related to PhoP), this protocol was initially used for vir- 
tual screening. 



Overview of experimental strategy 

Figure 3B illustrates an overview of the experimental strategy to 
assess biological activity and characterize the mode of action of 
the predicted compounds. As a primary assay, electrophoretic mobil- 
ity-shift assays (EMSA) were used to assess whether the com- 
pounds [n = 255) can inhibit the formation of the S. enterica PhoP- 
DNA complex necessary for gene regulation. To test whether the 
compounds active in EMSA [n = 8) affected PhoP dimerization, sec- 
ondary assays used included native polyacrylamide gel electrophore- 
sis (PAGE), size-exclusion chromatography, and Forster resonance 
energy transfer (FRET)-based assay. Native PAGE was used to 
detect qualitative effects on S. enterica PhoP dimerization. SEC was 
used to better characterize the compounds' ability to inhibit S. ent- 
erica PhoP dimerization by measuring the monomer-dimer profile. 
FRET was used as a higher-throughput method to quantify the 
effects of dimerization on £ coli OmpR/PhoB members to assess 
compound selectivity at the otA-fib-ab interface. Experimental results 
from the primary and secondary assays suggest the compounds 
inhibit formation of the PhoP-DNA complex, not by dimer inhibition, 
but in an allosteric manner to prevent DNA binding or by blocking 
the C-terminal DNA-binding domain. Potential mode of actions will 
be presented in a greater detail in Results and Discussions. 



Electrophoretic mobility-shift assays 

As a primary assay, an electrophoretic mobility-shift assay, also 
known as a gel-shift assay, was used to test the ability of the com- 
pounds to inhibit formation of the PhoP dimer-DNA complex neces- 
sary for gene expression. In the gel-shift assays, the protein is first 
activated by phosphorylation to induce the monomer-dimer equilib- 
rium. Radiolabeled DNA is incubated with the protein, and the com- 
pound of interest is subsequently added to the mixture. After an 
incubation period, the sample is loaded and electrophoresed on a 



gel to observe formation of the protein-DNA complex. Compounds 
that inhibit protein-DNA complex formation will lead to a darker 
band further down in the gel representing unbound DNA, while 
compounds that do not effect complex formation will result in a 
single band higher up on the gel (bound complex). Band intensities 
can be quantified and fitted to a sigmoidal curve to estimate IC 50 
values in dose-response experiments. 

PhoP was activated by in vitro phosphorylation with acetyl phosphate 
to induce the formation of monomer-dimer equilibrium. Compounds 
dissolved in DMSO (99.9%) were then incubated with the monomer- 
dimer mix, which did not exceed 1% volume of DMSO in the final 
solution concentration. Radioactive-labeled DNA containing phoP 
boxes (GGTTTAxxxxTGTTTA) were subsequently incubated with the 
mix to allow formation of the PhoP-DNA complex. Samples were 
loaded on 4-20% Tris-borate-EDTA (TBE) gels (Invitrogen) and electro- 
phoresed. Gels were dried and analyzed using a phosphorimager (Fuji- 
Film BAS-5000). Compounds that did not affect DNA binding 
displayed a single band upstream that represented the bound PhoP- 
DNA complex. Compounds that inhibited formation of the PhoP-DNA 
complex displayed a 'shift,' resulting in a band of the PhoP-DNA com- 
plex and/or a downstream band that represented unbound DNA. 

ImageJ 0 was used to quantitate the EMSA band intensities (for 
both bands representing bound and unbound DNA). KaleidaGraph^ 
was used to perform the curve fitting using a sigmoidal function to 
derive the IC 50 values from dose-response experiments. 

The p/jop-DNA fragments for electrophorectic mobility-shift assays 
were amplified by PCR using primers 312 and 369, and genomic 
DNA of wild-type S. enterica as template. The DNA fragments were 
isolated by running an electrophoretic gel and purified using QIA- 
quick columns (Qiagen). To radiolabel the DNA fragments, 100 ng of 
phop-DNA was used with T4 polynucleotide kinase and y- 32 P ATP 
and incubated at 37 °C for 1.5 h. Unincorporated DNA was 
removed using G-50 microcolumns (Amersham). 20 000 CPM of 
labeled probe (~12 fmol), 200 ng poly(dl-dC)-poly(dl-dC) (Amer- 
sham), and phosphorylated PhoP-His 6 were mixed with binding buf- 
fer (50 mini Tris-HCI pH 8.0, 50 mM KCI, 50 /.ig/mL BSA) in a total 
volume of 20 fil and incubated for 20 min at room temperature. 

For electrophoretic mobility-shift assays, S. enterica PhoP-His 6 was 
phosphorylated with acetyl phosphate. PhoP (0.6-1.2 nmol) was 
incubated in 50 (.il of phosphorylation buffer (50 mM Tris-HCI pH 
7.5, 50 mM KCI, 10 mM MgCI 2 ) containing 10 mM acetyl phosphate 
(Sigma-Aldrich) for 2.5 h at room temperature. Excess acetyl phos- 
phate was removed from phosphorylated PhoP-His 6 using a Micro 
Bio-Spin 6 Chromatography Column (Bio-Rad, Hercules, CA, USA) 
equilibrated with Tris buffer. Phosphorylated PhoP-His 6 was kept at 
4 °C and used within 24 h. 



Similarity search 

To identify similar compounds of the experimentally confirmed 
inhibitors, a similarity search (Tanimoto index) using the NCI website 
(http://129.43.27.140/ncidb2/) was performed to search the larger NCI 
library of 140 000 compounds. Using the experimentally confirmed 
compounds as the query, 36 compounds (Supporting Information) from 



Chem Biol Drug Des 2012; 79: 1007-1017 



1011 



Tang et al. 



the full NCI library were found, ordered, and experimentally tested by 
EMSA for inhibition activity. 



Native PAGE 

To test the effect of compounds on PhoP dimerization, a native 
polyacrylamide gel electrophoresis (PAGE) assay was used. The pro- 
tein is first activated by phosphorylation to induce formation of the 
monomer-dimer equilibrium (same method as described in EMSA). 
The compound of interest is incubated with the equilibrium mixture. 
Samples are loaded and electrophoresed on a native PAGE gel to 
separate the monomer and dimer states. Proteins can be visualized 
by gel staining. Compounds with no effects on dimerization should 
display two bands: one less intense and more downstream band 
representing the monomer, and a darker and more upstream band 
representing the dimer. Compounds with inhibition effects on dimer- 
ization will display a more intense band representing the monomer 
and a less intense band representing the dimer. 

Salmonella enterica PhoP-His 6 was phosphorylated in the same 
manner as by EMSA. Compounds were incubated with the active 
protein for 30 min, and loading buffer was added after. Samples 
were loaded and ran on a native gel (6% Tris-glycine; Invitrogen, 
Grand Island, NY, USA) at 4 °C. Gels were then stained with coo- 
massie blue to visualize the bands representing the protein. 



Size-exclusion chromatography (SEC) 

Size-exclusion chromatography was used to better characterize the 
monomer-dimer distribution. Protein is activated by phosphorylation 
(similar manner to EMSA and native PAGE) to induce the mono- 
mer-dimer equilibrium. Compound of interest is incubated with the 
equilibrium mix and injected to and separated by the size-exclusion 
column. Relative absorbance (A 2 8o) is detected from the elutions to 
quantify the amount of protein present. Compounds that do not 
inhibit dimerization will display a larger peak in the earlier elution 
representing the dimer, and a smaller peak in a later elution repre- 
senting the monomer. Compounds that do not effect dimerization 
will not display any changes in the monomer-dimer absorbance 
profile. Purified PhoP-His 6 was phosphorylated using 50 itim 
ammonium phosphoramidate [synthesized by the method of 
Sheridan ef al. (42)] and 20 mM MgCI 2 for 30 min at room tempera- 
ture. Compounds were subsequently incubated with PhoP for 
30 min. Samples of inactivated and activated PhoP (100 fil) were 
individually applied to a Superdex 75 column (GE Healthcare, Piscat- 
away, NJ, USA) equilibrated with elution buffer (50 mM Tris/HCI, 
pH 8.0, 150 mM KCI). Proteins were eluted in the same buffer at a 
flow rate of 0.5 mL/min. Protein concentration was assessed by 
measuring the 0D 280 . Fractions were collected and analyzed using 
native PAGE. A molecular weight standards kit (Sigma-Aldrich, St. 
Louis, MO, USA) was applied to the column for calibration. 



Forster resonance energy transfer analyses 

Forster resonance energy transfer (FRET) analyses were used to 
detect and quantify effects on dimerization on £ coll PhoP. In FRET, 
a cyan fluorescent protein-fused PhoP (CFP-PhoP) and a yellow fluo- 
rescent protein-fused PhoP (YFP-PhoP) are activated by phosphoryla- 



tion. Phosphorylation induces dimerization between CFP-PhoP and 
YFP-PhoP and brings the CFP and YFP in proximity where an energy 
exchange occurs. PhoP dimerization is characterized by a decrease 
in cyan emission and an increase in yellow emission. The rate of 
FRET increase depends on the rates of phosphorylation and dimer- 
ization. Compounds that inhibit dimerization will lead to a decrease 
in the FRET signal, while compounds that do not effect dimerization 
will not change the FRET signal. 

Forster resonance energy transfer analyses of fluorescent protein 
(FP) fused PhoP proteins, CFP-PhoP and YFP-PhoP, were performed 
as described in Gao ef al. (43). Phosphorylation of FP-PhoP was ini- 
tiated by addition of MgS0 4 and phosphoramidate to the mixture of 
CFP-PhoP, YFP-PhoP, and indicated compounds. The final concentra- 
tions were 0.6 fM CFP-PhoP, 2.5 (M YFP-PhoP, 20 mM phosphorami- 
date, 5 mM MgS0 4 , 100 fM compounds, and 1% (v/v) DMS0. 
Fluorescence was followed at 475 and 530 nm with excitation at 
430 nm. The ratio of 475 and 530 nm emissions was defined as 
the FRET ratio to monitor the interaction between CFP-PhoP and 
YFP-PhoP. 



Salmonella enterica PhoP expression and 
purification 

Salmonella enterica PhoP response regulator with a C-terminal His B 
tag was overexpressed in £. coli strain BL21 (DE3) transformed with 
the pT-7-7 plasmid. Cells were grown in Luria-Bertani medium with 
ampicillin (100 mg/L) and incubated until mid-logarithmic phase at 
37 °C. Overexpression was induced with 1 mM isopropyl-/J-D-thioga- 
lactopyranoside (IPTG) and incubated overnight at 25 °C. 

For purification, cells were harvested by centrifugation and then 
washed and resuspended in PBS, and stored overnight in -80 °C. 
Cells were then thawed at 4 °C, suspended in lysis buffer (50 mM 
l\laH 2 P0 4 , pH 8.0, 300 mM NaCI, 10 mM imidazole) and lysed by 
sonication. Cells were then centrifuged and the cell lysate was col- 
lected as the supernatant. Cell lysate was applied to a Ni 2+ column 
(GE Healthcare). Unbound proteins were removed with elution with 
Buffer A (50 mM Tris-HCI pH 8.0, 500 mM NaCI, 10 mM imidazole). 
Bound proteins were eluted with a 0-100% gradient of Buffer A to 
Buffer B (50 mM Tris-HCI pH 8.0, 100 mM NaCI, 500 mM imidazole). 
Fractions were collected and analyzed by SDS/PAGE. Proteins were 
concentrated using Amicon Ultra-15 filters (Millipore) and stored in 
storage buffer (20 mM Tris pH 7.8, 100 mM KCI, 20% glycerol) at 
-80 °C. Protein concentrations were determined by measuring the 
0D 280 (NanoDrop spectrophotometer). Purity was assessed by 
SDS/PAGE. 



Results and Discussion 

Computational strategy identified inhibitors 
that potentially bind at the PhoP ct4-p5-x5 
interface 

A total of 259 compounds were tested experimentally. From the 
results obtained by AutoDock and the consensus scoring scheme, 
the following number of compounds include 119 compounds that 
used a search grid centered at Arg 111 and 60 compounds that 
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used a search grid centered at Arg 118. From the results obtained 
by docking using FRED (search area centered at Arg 111), 40 com- 
pounds were tested. 



Eight compounds inhibited formation of the 
PhoP-DNA complex 

A total of 255 compounds were tested by EMSA: 119 compounds 
from the Arg 111 binding site screen using AutoDock, 40 com- 
pounds from the Arg 111 binding site screen using FRED, 60 com- 
pounds from the Arg 118 binding site screen, and 40 compounds 
identified from a Tanimoto similiarity search for analogs of the ini- 
tial hits within the entire NCI library. Eight compounds (NSC9B08, 
NSC45576, NSC48B30, NSC35489, NSC65238, NSC88915, 
NSC1 18806, and NSC168197) (Figure 5) displayed inhibition activity 
via disruption of the PhoP-DNA complex. The eight compounds 
inhibited PhoP-DNA complex formation in a dose-dependent manner 
(Figure 6), with six of the eight compounds displaying an IC 5D in the 
micromolar range (from 3.6 to 285 /Ml). 

Analogs (36 in total) of the eight inhibitors based on Tanimoto simi- 
larity (90% similar) from the NCI library (140 000 compounds) did 
not display inhibition activity by EMSA. Inactivity of these analogs 
may be due to the sensitive preference of the PhoP binding site at 
the a4-/?5-a5 interface. 



Compounds did not effect PhoP dimerization 

The eight compounds did not display any effects on band intensity 
for the band representing the PhoP monomer and dimer in the 
native PAGE assay (Figure 7). These results suggested that the com- 
pounds do not act as dimerization inhibitors. However, band inten- 
sity changes may not have been observed because of the inability 
of the gel to distinctly separate the phosphorylated monomer and 
dimer species. 



The eight compounds also did not display any changes on the 
monomer-dimer profile of PhoP as compared with a control (DMSO) 
using size-exclusion chromatography, which performs a better 
separation and characterization of the oligomeric states than native 
PAGE. One caveat of using SEC as a binding assay is that a 
compound exhibiting a rapid off-rate cannot be detected with SEC 
because of time (30 min.) needed to perform the separations. 

Results from the FRET assays using £ coli PhoP also did not 
suggest any degree of dimer inhibition (Figure 8). However, with 
the use of CFP-PhoP and YFP-PhoP, one possibility is that the 
compounds may bind to CFP or YFP instead of PhoP, resulting in 
undetected inhibition effects. Some compounds (e.g., 1, 2) were 
inherently colored and affected the detection of FRET signal and, 
therefore, could not be properly tested for dimer inhibition. Com- 
pounds 1, 3, and 8 displayed a smaller FRET ratio changes as com- 
pared with the positive control sample (e.g., CY +DMS0, CY PhoP). 
However, compounds 1 and 3 interfere with the FP fluorescence, 
while compound 8 alters the fluorescence ratio of FP-PhoP pairs. 
Therefore, the FRET method is not sufficient to assess the inhibition 
of compounds 1, 3, and 8. FRET also cannot detect compounds that 
bind to the ocA-fib-ab interface to inhibit DNA binding in an 
allosteric manner, because only the N-terminal of the PhoP 
response regulator bound to CFP or YFP was used. 



Compounds suggested selective inhibition for 
Salmonella enterica PhoP 

The 8 compounds did not display any effects on the signal in the 
FRET assays with £ coli PhoP, which suggested no inhibition of 
dimerization (Figure 8). The compounds may bind to the £ coli PhoP 
a4-/?5-a5 interface, but have no effects on dimerization. However, 
because the £ coli PhoP linker and C-terminal were not included in 
the FRET analyses, there remains a possibility that the compounds 
may bind to £ coli PhoP to inhibit DNA binding. These results 
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Figure 5: Chemical structures of the eight compounds (1-8) that inhibited formation of the PhoP-DNA complex by electrophoretic mobility- 
shift assays (EMSA) with their NSC number and their estimated IC50 values. 
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Figure 6: Dose response curves (fitted by sigmoidal function, 
KaleidaGraph) from electrophoretic mobility-shift assays in presence 
of compounds 1-4 (A) and 5-8 (B). The points plotted represent 
the relative amount of PhoP-DNA complex formed in EMSA 
experiments. 

obtained so far suggested that the compounds might selectively 
bind to the S. enterica PhoP a4-/?5-a5 dimer interface to inhibit 
formation of the protein-DNA complex, perhaps by an allosteric 
mechanism. 



Inhibitors may act by a different mode of 
action 

Experimental results from a series of in vitro and biophysical assays 
(EMSA, Native PAGE, SEC, FRET) suggested that the compounds 
tested disrupt PhoP-DNA complex formation, but not via direct ho- 
modimer inhibition. Instead, one potential mode of action is to bind 
to the C-terminal DNA-binding domain, a more direct mechanism to 
disrupt formation of the PhoP-DNA complex. Another possible mode 
of action is by binding to the N-terminal regulatory receiver domain 
or the linker region to act in an allosteric manner and prevent con- 
formational changes necessary for DNA binding. The former may 
prove to be a more effective strategy for selective inhibition of 
response regulators, because of the conserved aA-pb-xb structural 
motif of the dimer interface. 

Virtual screening may have helped to identify compounds binding 
the interface. The inhibitors did not block dimerization, but did pre- 
vent formation of PhoP-DNA complex necessary for gene regulation, 
perhaps by binding at the cx4-/?5-a5 interface and acting by an allo- 
steric mechanism. 



Challenges of structure-based design to target 
protein-protein interactions 

In the pursuit of discovering PPI inhibitors, the original objective 
was to use small molecules as probes to identify and character- 
ize potential binding sites at the a4-/?5-a5 interface of S. enteri- 
ca PhoP to block homodimerization. Through the discovery of 
dimerization inhibitors presumably binding at the sites where the 
critical residues are located, it was hoped that the findings 
would lead to a better understanding of the physical properties 
underlying molecular recognition of protein-ligand interactions at 
PPI. While the original intent of this study was to target the 
o£4-/?5-o:5 interface to block dimerization, the experimental results 
suggest an unexpected finding: drug-like compounds may bind at 
the cx4-/85-a5 interface and function in an allosteric manner and 
cause a conformational change to prevent DNA binding. Further 
experimental characterization (e.g., crystallography, NMR) will be 
necessary to identify potential binding sites and to elucidate the 
atomic details of the complexes to determine the mode of 
action. 

Structure-based discovery of protein-protein interaction inhibitors 
remains a significant challenge as seen by the high percentage of 
false positives from the computational predictions in this study. 
Docking and structure-based design methods to incorporate protein 
flexibility will need to be used to identify and design potential 



Acetyl Phosphate - + + + + + + f + + + 
Compound (1 00 tm) - - DMSO 1 2 3 4 5 6 7 8 

Figure 7: Native PAGE results of S. enterica PhoP in presence of eight inhibitor compounds (1-8). PhoP dimerization was induced by 
phosphorylation via acetyl phosphate. An unphosphorylated sample without compound (first lane from left) was used as a negative control. 
Samples with buffer only (second lane from left) and with 1% (v/v) DMSO (third lane from left) were used as positive controls. 
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Figure 8: Time-dependent FRET of Escherichia coli FP-PhoP pairs 
in the presence of inhibitor compounds (1-8). Inactive PhoP (CY- 
PhoP no P) was used as a negative control, while activated PhoP 
(CY-PhoP) and sample with DMSO (CY +DMS0) were used as posi- 
tive controls. 

allosteric inhibitors or ones that bind by an induced-fit mechanism. 
Another potential limitation is in the scoring functions used for 
virtual screening, which cannot accurately predict binding affinities, 
especially for systems that are not present in the scoring function 
training sets (e.g., binding sites with relatively flat surfaces). Even 
if the docking modes may be correctly identified, limitations in scor- 
ing accuracy may classify potential tight-binding ligands as weak- 
binders, resulting in a high percentage of false positives. Perhaps, 
the use of first-principle methods for estimating binding affinities 
may help circumvent this limitation, assuming that the correct 
ligand conformation was predicted by the docking procedure. 

Also, availability of a high-resolution crystal structure of PhoP 
bound to its DNA promoter should better elucidate the conforma- 
tion of the «4-j85-a5 interface in the biologically active form bound 
to DNA and provide atomic-level detail of potential binding sites 
necessary for docking and structure-based drug design. However, 
inherent inaccuracy of using a 'static' structure for molecular design 
may lead to false negatives using structure-based design methods 
without incorporation of protein flexibility. For more accurate com- 
putational modeling, in particular to target PPI, improvements in 
both docking and scoring methods are necessary. With the increas- 
ing interest in targeting PPI and availability of structural and bind- 
ing affinity data, development of more accurate and robust SBDD 
methods to target PPI will become possible. 



Discovery of first-in-class two-component 
signal transduction response regulator 
inhibitors 

In this study, eight first-in-class inhibitors of the S. enterica PhoP 
TCST response regulator were discovered using a hybrid approach 



coupling computational and experimental methods for molecular 
design. Potential mode of action of the compounds was character- 
ized by a series of in vitro and biophysical assays. Compounds may 
potentially bind at the a4-/i5-a5 interface and act as allosteric 
inhibitors, rather than dimerization inhibitors, to prevent DNA bind- 
ing. Discovery of first-in-class PhoP inhibitors should serve as a 
proof-of-concept for targeting TCST response regulators as a novel 
strategy to inhibit bacterial virulence. 



Conclusions and Future Directions 

Targeting two-component signal transduction response regulators 
to modulate virulence gene expression is a promising strategy 
for antibiotics development. With the increasing resistance of 
bacterial pathogens to current therapeutics, antibiotics that 
can prevent virulence instead of inhibiting growth or inducing 
death may lead to less selective pressure for the generation of 
resistance. 

In this study, eight compounds have been discovered by coupling 
computational and experimental methods to disrupt formation of 
the S. enterica PhoP-DNA complex necessary for gene regulation. 
Eight compounds inhibited the PhoP-DNA complex formation in a 
dose-dependent manner by EMSA. Based on the experimental 
results, the PhoP inhibitors may potentially bind to the plastic 
oc4-/J5-ot5 interface and act by an allosteric mechanism to 
prevent DNA binding. Alternative modes of action include binding 
to other regions of the N-terminal domain to act in an allosteric 
manner or the C-terminal DNA-binding domain to directly inhibit 
formation of the PhoP-DNA complex. Experimental results 
obtained from a series of in vitro and biophysical assays suggest 
a potential of these compounds to inhibit bacterial virulence. 

Further elucidation of the mode of action to assess the potential of 
the eight compounds as virulence inhibitors is planned. Structural 
analogs of the eight compounds can be designed to enhance affin- 
ity and characterize structure-activity relationships. Structural stud- 
ies such as X-ray crystallography and NMR must be performed to 
validate the mode of action and elucidate the protein-ligand inter- 
actions in atomic detail. 
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